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ABSTRACT 

Sensitive 5-38 /im Spitzer Space Telescope (SST) and ground based 3-5 /zm spectra of the embedded 
low mass protostars B5 IRS1 and HH46 IRS show deep ice absorption bands superposed on steeply 
rising mid-infrared continua. The ices likely originate in the circumstellar envelopes. The CO2 bending 
mode at 15 fiin is a particularly powerful tracer of the ice composition and processing history. Toward 
these protostars, this band shows little evidence for thermal processing at temperatures above 50 K. 
Signatures of lower temperature processing are present in the CO and OCN~ bands, however. The 
observed CO2 profile indicates an intimate mixture with H2O, but not necessarily with CH3OH, in 
contrast to some high mass protostars. This is consistent with the low CH3OH abundance derived 
from the ground based L band spectra. The CO2/H2O column density ratios are high in both B5 
IRS1 and HH46 IRS (~ 35%). Clearly, the SST spectra are essential to study ice evolution in low 
mass protostellar environments, and to eventually determine the relation between interstellar and 
solar system ices. 

Subject headings: Infrared: ISM — ISM: molecules — ISM: abundances — stars: formation — stars: indi- 
vidual (B5 IRS1) — astrochemistry 



1. INTRODUCTION 

A recurring question in disk, planet and comet forma- 
tion studies is how the composition of molecular material 
evolves as it flows from the molecular cloud to the pro- 
tostellar envelope and protoplanetary disk. Much of the 
material in these environments is frozen on grains. A 
plethora of processes, including heating by stellar pho- 
tons, shocks related to accretion or outflow, cosmic ray 

*SOME OF THE DATA PRESENTED HEREIN WERE OBTAINED AT THE 
W.M. KECK OBSERVATORY, WHICH IS OPERATED AS A SCIENTIFIC PART- 
NERSHIP AMONG THE CALIFORNIA INSTITUTE OF TECHNOLOGY, THE 
UNIVERSITY OF CALIFORNIA AND THE NATIONAL AERONAUTICS AND 
SPACE ADMINISTRATION. THE OBSERVATORY WAS MADE POSSIBLE BY 
THE GENEROUS FINANCIAL SUPPORT OF THE W.M. KECK FOUNDATION. 

'THE VLT/ISAAC SPECTRA WERE OBTAINED AT THE EUROPEAN SOUTH- 
ERN OBSERVATORY, PARANAL, CHILE, WITHIN THE OBSERVING PRO- 
GRAM 272.C-5008 

1 Division of PMA, Mail Code 105-24, California Institute of 
Technology, Pasadena, CA 91125, USA 

2 Leiden Observatory, PO Box 9513, 2300 RA Leiden, the 
Netherlands 

3 SRON, PO Box 800, 9700 AV Groningen, the Netherlands 

4 Division of GPS, Mail Code 150-21, California Institute of 
Technology, Pasadena, CA 91125, USA 

5 Max-Planck-Institut fur Astrophysik, P.O. Box 1317, D-85741 
Garching, Germany 

6 Department of Astronomy, University of Texas at Austin, 1 
University Station C1400, Austin, TX 78712-0259 

7 Spitzer Science Center, California Institute of Technology, CA 
91125, USA 

8 Smithsonian Astrophysical Observatory, 60 Garden Street, 
MS42, Cambridge, MA 02138 

9 Department of Physics and Astronomy, Northern Arizona Uni- 
versity, Box 6010, Flagstaff, AZ 86011-6010 

10 Department of Astronomy, University of Maryland, College 
Park, MD 20742 

Jet Propulsion Laboratory, MS 183-900, California Institute 
of Technology, 4800 Oak Grove Drive, Pasadena, CA 91109 



hits, and ultraviolet irradiation can change the ice struc- 
ture and composition. The spectroscopic effects of these 
processes can be observed in the vibrational bands of 
the ices through infrared spectroscopy. Ices around high 
mass protostars have been extensively studied this way. 
Infrared Space Observatory (ISO) spectra have shown 
that in particular the ice structure is affected by heating 
from the central star. The simplicity of the ice compo- 
sition does indicate that the formation of new species 
through ultraviolet irradiation or cosmic ray hits occurs 
at a low level at best. Observations of ices around low 
mass protostars have been limited due to the unavailabil- 
ity of much of the 5-20 /im spectral region, where many of 
the molecular bending mode transitions occur. In par- 
ticular the CO2 bending mode at 15 /xm is a valuable 
tracer of ice structure and composition (Ehrenfreund et 
al. 1998; Gerakines et al. 1999). With the se nsitive In- 
frared Spectrometer (IRS: iHouck et alj|2004l) on board 
of the Spitzer Space Telescope fSST: lWerner et al.ll2004|) 
this band can now be observed for the first time at high 
quality in low mass systems. 

Observation of two protostars are presented in this pa- 
per. B5 IRS1 (IRAS 03445+3242; Beichman et al. 1984; 
L =10 Lq) is well studied at infrared and millimeter 
wavelengths (e.g. Charnley et al. 1990; Langer et al. 
1996). The millimeter continuum emission is resolved on 
a few arcsec scale, and may originate in an inclined disk. 
The outflow has received most of the attention, and has 
a large opening angle leading to significant outflow/infall 
interaction. HH46 IRS (IRAS 08242-5050; L = 12 L Q ) 
is also deeply embedded, and is also the driving source 
of a powerful outflow. SST imaging and spectroscopic 
observations of this source, focused on the outflow, are 
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Fig. 1. — Combined Spitzcr Space Telescope and ground based L 
and M band spectroscopy of B5 IRS1 (top, multiplied by factor of 
5 for clarity) and HH46 IRS (bottom). Identifications and possible 
identifications are indicated. 



presented in Noriega-Crespo et al. (2004). 
2. OBSERVATIONS 

B5 IRS1 and HH46 IRS were observed with SST/IRS 
as part of the 'c2d' Spitzer Legacy program (Evans et 
al. 2003) in the modules Short Low (SL; A = 5 — 14 /xm; 
R=64-128), Short High (SH; A = 10-20 /mi; R=600) and 
Long High (LH; A = 20 - 38 /im; R=600). The archival 
AOR keys are 0005638912 (HH46 IRS) and 0005635328 
(B5 IRS1) for PROGID 172. Both sources are well cen- 
tered in all slits. The integration time was 28 seconds 
per module per source at 14 second ramps, except for 
SH which has 24 seconds in total and 6 second ramps. 
The spectra were reduced with the IRS pipeline version 
9.5 on 10 March 2004 at the SSC. Bad pixels were in- 
terpolated in the spectral domain in the 2D images be- 
fore extracting ID spectra. Accurate wavelength cal- 
ibration was assured using calibration tables available 
in May 2004. For overlapping SH spectral orders, the 
poorly calibrated long wavelength part of each order was 
removed. This is particularly important to obtain a re- 
liable profile of the CO2 ice band at 15.2 /im, where two 
orders overlap. HH46 IRS was observed independently 
as an Early Release Observation (AOR key 0007130112 
and PROGID 1063; Noriega-Crespo et al. 2004). The 
datasets are in good agreement, and were averaged. Fi- 
nally, complementary ground based observations were 
obtained. B5 IRS1 w as observed with Keck/NIRSPEC 
IjMcLean et al. 111998(1 at R=25000 in t he M band and at 
R=20 00 in the L band. VLT/ISAAC (jMoorwood et all 
1999) L and M band spectra of HH46 IRS were obtained 
at R=600 and 5000 respectively. The SST and ground 
based spectra were combined by scaling the ground based 
data. 

3. RESULTS 

The combined SST and ground based 3-38 ^m spectra 
of B5 IRS1 and HH46 IRS show steeply rising continua 
(Fig. 1). The continuum of HH46 IRS is the steep- 



est, with a 35 /xm/4 /xm flux ratio of 100, compared to 
10 for B5 IRS1. Numerous silicate and ice absorption 
features are superposed. The main ice constituents are 
H2O, CO2, and CO. The well known, yet unidentified, 
3.47 and 6.85 /xm bands are present in both sources as 
well. Rarely seen before is a shallow absorption feature 
at 11.2 /im in both sources. It may be related to re- 
fractory dust compo nents, such as crystalline forsterite 
l|Kessler et al.ll2004D . In addition, HH46 IRS, but not 
B5 IRS1, shows absorption by CH 4 and 'XCN' (likely 
OCN - ; van Brockhuizen et al. 2004). Finally, other ice 
features may be present in the 5-10 /xm region (NH3, 
HCOOH), but a dedicated analysis is required to verify 
their reality. We focus on the band profiles and abun- 
dances of the main ice components and their relation to 
the 15 /im CO2 bands, newly discovered with the SST. 

The CO2 bending mode has rarely been observed to- 
ward low mass protostars, and never at such high quality. 
It is blended with the short wavelength wing of the sili- 
cate bending mode. The CO2 band is put on an optical 
depth scale, assuming the 'intrinsic' profile of the silicate 
band is represented by the Galactic Center source GC3 
(Chiar et al. 2000). For this a third order polynomial 
is fitted to the wavelength regions 13.0-14.7 and 26.3- 
33.3 /im. The resulting CO2 bands look similar to those 
observed in other lines of sight; note in particular the 
presence of a long wavelength wing extending to at least 
16 /xm (Gerakines et al. 1999). Similar to the massive 
protostar NGC 7538 IRS9 (Fig. 2b), weak double ice 
crystallization peaks are observed at the bottom of the 
band in HH46 IRS, but not in B5 IRS1. Neither B5 IRS1 
nor HH46 IRS show evidence for a third peak at 15.38 
/xm, expected in CH30H:C02 complex formation. 

To place the interpretation of the CO2 band within 
a larger perspective, ground based observations of the 
3.53 /xm CH3OH band and the 4.67 /xm band of solid 
CO are analyzed. The 3.53 /xm band is superposed on 
the wing of the strong 3.07 /xm H2O band and is locally 
blended with the unidentified 3.47 /xm band. Following 
Brooke et al. (1999), we derive the continuum and sep- 
arate the CH3OH contribution. For comparison, spectra 
of the massive protostars W33A (Brooke et al. 1999) and 
NGC 7538 IRS9 (Boogert et al. 2004) are analyzed as 
well. Hints of CH 3 OH are seen in both HH46 IRS and B5 
IRS1, resulting in column densities of 7% with respect to 
H2O, comparable to NGC 7538 IRS9, but much less com- 
pared to W33A. For HH46 IRS the detection of CH 3 OH 
is strengthened by the presence of a feature at 9.7 /xm 
in the bottom of the silicate band (Fig. 1). Both HH46 
IRS and B5 IRS1 show prominent bands of solid CO at 
4.67 /xm (Fig. 2d). The ratio between the central nar- 
row CO component and the broad long wavelength wing, 
representing the column density ratio of volatile pure CO 
and CO mixed with less volatile H2O, is a factor of five 
smaller in HH46 IRS. In fact, the profile of HH46 IRS is 
comparable to that of the massive protostar W33A. The 
latter two sources also show a band at 4.62 /xm, most 
likely attributed to the OCN - ion. Relative to H2O, the 
OCN - column density is comparable to (upper limits to) 
those of B5 IRS1 and NGC7538 IRS9. Furthermore, the 
4.67 /xm spectrum of B5 IRS1 shows deep ro- vibrational 
gas phase CO lines. The presence of gas phase CO in 
HH46 IRS is hard to assess because of telluric contamina- 
tion. The intriguing differences and similarities between 
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Fig. 2.— Solid CH 3 OH (a), CO2 (b and c), and CO (d) spectra of HH46 IRS and B5 IRS1 compared to the massive protostars NGC7538 
IRS9 and W33A in each panel from top to bottom. The normalization factor is indicated in the right corner of each plot. Panel a: a 
laboratory pure CH3OH spectrum at T lab = 10 K (thick yellow line) to which is added a Gaussian (dashes) to account for the underlying 
"3.47 /urn" absorption feature. The NGC7538 IRS9 and W33A data are from Boogert et al. (2004) and Brooke et al. (1999) respectively. 
Panel b: the C0 2 bending mode observed with SST/IRS (HH46 IRS and B5 IRS1) and ISO/SWS (NGC 7538 IRS9 and W33A). The 
thick yellow dashed line is the laboratory ice H2O:CH3OH:CO2=0. 84:0. 92:1, and the continuous thick yellow line is a laboratory ice with 
reduced CH3OH content, H 2 O:CH 3 OH:CO2=0.92:0.29:l, both at T lab =115 K. Only for W 33A an enhanced CH3OH/CO2 ratio is required, 
in agreement with the strong 3.53 /im CH3OH feature in panel (a). Panel c: the CO2 bending mode for the same sources as in panel 
(c), but now fitted with different laboratory spectra: a highly processed polar ice (Ti ab =125 K; red) in combination with an H20-rich, 
CHsOH-dcficicnt cold ice (Ti ab =10 K; blue). The sum of these two components is the thick yellow line. Panel d: the absorption band 
of solid CO observed with VLT/ISAAC (HH46 IRS), Keck/NIRSPEC (B5 IRS1 and NGC7538 IRS9) and ISO/SWS (W33A). Narrow 
absorption lines are due to circumstellar gas phase CO. The onset of the 4.62 fim 'XCN' band is seen in HH46 IRS and W33A. Gaussian 
and Lorentzian fits to the apolar and polar CO components are indicated by blue and red lines respectively. The thick yellow line represents 
the combined fit, including the 'XCN' bands and blue-shifted CO component that are not shown separately. 



HH46 IRS and B5 IRS1, as well as compared to mas- 
sive protostars, provide insight into the formation and 
evolution of interstellar and circum-protostellar ices. 

Finally, column densities of the main ices, summarized 
in Table 1, are derived by dividing the integrated optical 
depth over the laboratory integrated band strength (e.g. 
Hudgins et al. 1993). Note that the C0 2 /H 2 ratios 
toward B5 IRS1 and HH46 IRS are significantly larger 
compared to the average over many, mostly massive pro- 
tostellar sightlines (0.17±0.03; Gerakines et al. 1999). 

4. DISCUSSION 



4.1. Evolution of Ices in Low Mass Environments 

The formation and evolution of interstellar ices is, in 
principle, strongly dependent on local conditions such as 
the atomic hydrogen and carbon density, the tempera- 
ture, the cosmic ray flux, the ultraviolet photon flux, and 
the presence of shocks. Thus, key issues are the location 
of the ices along the absorption line of sight and the rela- 
tive contributions from foreground clouds, envelopes, and 
inclined disks. The continuous rise of both the B5 IRS1 
and HH46 IRS SEDs between 3 and 40 ^m, as well as the 
detection of extended submillimeter emission in JCMT 
archive images, are in favor of envelope-dominated mod- 
els. We established the properties of these envelopes us- 
ing the approach of J0rgensen et al. (2002). Adopting 
optical constants for bare and ice coated silicate grains 



(Ossenkopf & Henning 1994) the 2-2000 /im SED and the 
depth of the observed superposed ice and silicate bands 
were self-consistently modeled. The SED and silicate 
band depth of B5 IRS1 are well fitted by the envelope 
models, and a significant fraction of the ices has evap- 
orated. In contrast to B5 IRS1, the submillimeter/far- 
infrared and mid-infrared SST SED of HH46 IRS cannot 
be simultaneously fitted. Possibly, this envelope is em- 
bedded in a larger scale, cold cloud not modeled within 
the framework of the simple spherical envelope. 

Next we address the extent to which the ices in these 
envelopes have been processed. Several indicators are 
available. Laboratory experiments have shown that heat- 
ing of ice mixtures with concentrations C02/H20> 1 
results in crystallization and an effective segregation of 
the CO2 and H2O species. Spectroscopically this is rec- 
ognized as double peaked profiles, characteristic of the 
pure CO2 matrix (Ehrenfreund et al. 1998). Depend- 
ing on the ice composition, the amorphous to crystalline 
phase transition occurs at 50-90 K in space, lower than 
the corresponding laboratory temperatures due to the 
longer interstellar time scales (Boogert et al. 2000). Sub- 
structures are seen in the HH46 IRS CO2 band (Fig. 2b), 
but are much weaker compared to some highly heated 
massive protostellar envelopes (Gerakines et al. 1999), 
and are absent in B5 IRS1. In fact, the CO2 band 
profile of HH46 IRS is comparable to that observed to- 
ward one of the least processed envelopes, surrounding 
the massive protostar NGC 7538 IRS9. In the simplest 
scenario of a single ice at one temperature, the mixture 
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TABLE 1 

Column densities and Abundances 



Quantity 



B5 IRS1 HH46 N7538/9° W33A f 



N(C0 2 )/N(H 2 0) 

Af(CO-p)/Af(H 2 0) a 

A/"(CO-np)/AT(H 2 0) b 

Ar(CH 3 OH)/Af(H 2 0) 

7V(OCN-)/AT(H 2 0) c 



0.37 


0.32 


0.24 


0.16 


0.16 


0.15 


0.02 


0.08 


0.27 


0.05 


0.14 


0.03 


<0.06 


0.07 


0.07 


0.22 


< 0.005 


0.007 


0.006 


0.02 



AT(H 2 0) [10 18 cm" 2 ] 

N(H 2 0)/N H [10- 5 ] 
N H [10 22 cm" 2 ] 



2.2 (0.3) 8.0 
5.2 5.7 

4 d 14 d 



6.8 

5.0 
16 



9.0 

3.2 

28 



a broad 'polar' CO component, likely H 2 mixture 

b broad 'apolar' CO component, likely pure CO 

c using OCN - band strength of van Broekhuizcn et al. 2004 

d from the envelope models described in the text 

c Whittct et al. 1996 and references therein 

f Gibb et al. 2000 and references therein 



CH 3 OH:H 2 O:CO 2 =0.3:l:l (§4.2) must still have a labo- 
ratory temperature as high as 115 K, or ^75 K in space 
(Fig. 2b). In the more likely scenario of a temperature 
gradient along the line of sight, the bulk of the HH46 
IRS and B5 IRS1 envelopes have temperatures well be- 
low 50 K. A fraction of the inner envelope of HH46 IRS 
must be warmer, causing the observed weak substruc- 
tures. Such two component fits explain the observed 
profiles well (Fig. 2c). 

More extensive processing at lower temperatures most 
likely has occurred within the envelopes, however. Unlike 
the massive protostar NGC 7538 IRS9, the ground based 
4.67 /im spectra of solid CO toward HH46 IRS and B5 
IRS1 show a weak central narrow component (Fig. 2d; 
Table 1). The most volatile CO-rich 'apolar' ices may 
thus have evaporated (Tielens et al. 1991). HH46 IRS 
shows a particularly broad profile, resembling the mas- 
sive protostar W33A. Like W33A, HH46 IRS shows an 
absorption feature at 4.62 /im, likely attributed to the 
OCN - species. This molecule may be produced from 
HNCO in the solid state at relatively low ice tempera- 
tures <50 K (van Broekhuizen et al. 2004). 

Concluding, while high temperature ice processing, 
traced by the CO2 band, is not observed in the low 
mass envelopes, low temperature (< 50 K) processing 
may play a significant role. Qualitatively this is similar 
to some high mass protostars (W33 A, NGC7538 IRS9). 
The proposed evolutionary sequence of ice processing in 
massive envelopes (Boogert et al. 2000; van der Tak et 
al. 2000) is indeed largely based on high temperature in- 
dicators, such as ice crystallization, hot core gas temper- 
atures and gas/solid state ratios. These relations need to 



be investigated in a larger sample of low mass protostars, 
in which the CO2 band profile is a crucial tracer. 



4.2. CH 3 OH:C0 2 Complexes 

The long wavelength wing of the CO2 bending mode 
may be a tracer of the presence of CH3OH in the ices. 
CO2 and CH3OH form complexes, leading to an en- 
hanced wing in some sources well fitted by laboratory 
ices with a CH30H:H"20:C02 mixing ratio of 1:1:1 (e.g. 
W33A in Fig. 2b). Recently, large abundances of 
CH3OH of 25% with respect to H 2 were found in the 
envelopes of some low mass protostars (Pontoppidan et 
al. 2003). Clearly, B5 IRS1 and HH46 IRS both have 
lower CH3OH abundances (Table 1). This is consistent 
with the weakness of the long wavelength CO2 wing. 
Indeed, laboratory mixtures with a CH30H:H20:C02 
ratio of 0.3:1:1 fit the observed band well. The high 
CO2/CH3OH column density ratio of ^6 toward these 
sources is however barely consistent with this laboratory 
mixture. Alternatively, the band profile can be explained 
by the combination of an abundant H^O-rich ice respon- 
sible for the long wavelength wing, and an at least partly 
heated C02-rich ice responsible for the crystallization 
substructures seen in HH46 IRS (Fig. 2c). 



5. CONCLUSIONS AND FUTURE WORK 

High quality SST observations of the CO2 bending 
mode at 15 /jm toward low mass protostars offer a new 
tracer of ice mantle composition and evolution. The em- 
bedded low mass systems B5 IRS1 and HH46 IRS show 
CHsOH-poor ices with little evidence for 50-90 K ther- 
mal processing in their envelopes. Lower temperature 
processing appears evident in the solid CO band. These 
results form the basis for future studies on the physical 
and chemical state of ices entering protoplanetary disks, 
and how these and solar system ices are related. CO2 
bending mode observations of more evolved systems and 
edge-on disks are required. 
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